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The evolution of the red sequence slope in massive galaxy clusters 
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ABSTRACT 

We investigate the evolution of the optical and near-infrared colour-magnitude relation in an 
homogeneous sample of massive clusters from z = 1 to the present epoch. By comparing deep 
Hubble Space Telescope ACS imaging of X-ray selected MACS survey clusters at z ~ 0.5 to 
the similarly selected LARCS sample at z ~ 0.1 we find that the rest-frame S(U — V)/SV 
slope of the colour-magnitude relation evolves with redshift which we attribute to the build up 
of the red sequence over time. This rest frame slope evolution is not adequately reproduced 
by that predicted from semi-analytic models based on the Millennium Simulation despite a 
prescription for the build up of the red sequence by in-falling galaxies, 'strangulation' . We ob- 
serve no strong correlation between this slope and the cluster environment at a given redshift 
demonstrating that the observed evolution is not due to a secondary correlation. Also pre- 
sented are near-infrared UKIRT WF CAM obse r vation s of the LARCS clusters which confirm 
and improve on the the result from lStott et alj {2007) finding that there has been a two-fold 
increase in faint My > — 20 galaxies on the red sequence since z = 0.5 to a significance of 
5a. 

Key words: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - galaxies: 
evolution 



1 INTRODUCTION 

Galaxy clusters are important laboratories for the study of galaxy 
formation and evolution as they contain a concentrated population 
of diverse galaxies in a relatively small volume. Early workers in 
the field found that if a colour-magnitude diagram is plotted for 
members of local clusters such as Virgo and Coma, the Elliptical/SO 
galaxies wer e found to be confined to a prom inent linear feature in 
colour space dVisvanathan & SandagJ 19771) . This feature is known 
as the red sequence and has a very small intrinsic scatter (typically 
<0.1 mag) which has been interpreted as evidence that the passive 
galax ies in clusters formed coevally at high redshift (Bow er et al.1 
|l993). 

The red sequence is observed to have a slope such that the 
faint galaxies are bluer than the bright cluster members. The origin 
of the slope has been controv ersial due to the age m etallicity degen- 
eracy for stellar populations. IWorthev e t al. ( 1995) showed that the 
sequence of colours that comprise the slope can be equally well ex- 
plained by a progressive decrease in either metallicity or stellar age. 
The slope is now thought to be due to a mass-metallicity relation 
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along the red sequence as this degeneracy was broken by compar- 
ing co lour-magnitude simulatio ns to observations of distant clus- 
ters dKodama & Arimotoll 19971) . The origin of the mass-metallicity 
relation is thought to be the heating of the interstellar medium 
(ISM) by supernovae which triggers the formation of a galactic 
wind when the thermal energy of the gas exceeds the binding en- 
ergy. This wind ejects gas more efficiently in smaller galaxies due 
to their shallower potential wells, resulting in the trend of increased 
metallicity with mass which manifests itself as the massive metal- 
rich galaxies appearing progressively redder than their less massive 
counterparts (Carlberg Hl984lArimoto & Yoshiilll987t) . Within the 
hierarchical picture of galaxy evolution, red sequence galaxies form 
from the mergers of star forming disc galaxies with the largest el- 
lipticals being the product of the merger of the largest disc sy stems 
and are thus the most metal rich ( Kau ffmann & Charlotll 19981) . 

The observed frame colour-magnitude relations for clusters 
at similar redshifts are found to have compar able slopes which 
have been shown to change with redshif t dLopez-Cruzl 1 19971 : 
iGladders et alJll998l ; lLopez-Cruz et ai]|2004l) . One explanation for 
this change in the observed slope is that it is a result of K cor- 
rection and an evolution in the mass m etallicity relation in cluster 
ellipticals ( Koda ma & Arimotol 1 19971) . The change of slope with 
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redshift has been used as a method to constrain cluster galaxy evo- 
lution, with comparisons of the observed slope evolution to models 
suggesting that the elliptical galaxies in cluster co res have been in 
place since at least z = 2 dGladders et al .1 1998). However other 
studies of the rest frame slope have found little o r no evolution 
suggesting that K correction is the dominant factor dStanford et al.l 
1 19981 : iBlakeslee et all [2003; Meietal] 120061) . This non-evolution 
has been used as evidence to favour monolithic collapse over hier- 
archical methods of galaxy formation. 

Current research on the build up of the red sequence suggests 
there may indeed be an age contribution to the red sequence slope 
and its evolution as faint, sub-L* galaxies, thought to have under- 
gone recent star formation, are found to be transforming onto the 
sequence over time as they fall into the dense cluster env ironment 
dPe Lucia et ai]|2007l ; [stott et al.ll2007l ; ISmifh et al.ll2008h - This is 
thought to be connected to the observation that there is a greater 
abundance of SO galaxies in local clusters compared to their high 
redshift counterparts, whose progenitors may be the quenched star 
forming galaxies in higher redshift systems (Dressie r et aljfl997l : 
Smail et al. 1998). If this transformation scenario is correct then 
we expect to see evidence for it in the evolution of the rest frame 
colour-magnitude relation. 

Regardless of its origin, empirical observations and models 
of the colour and slope of the observed red sequence can be used 
in combination as a template to predict the redshift of a cluster. 
This is found to be in good agreement with spectroscopic red- 
shifts and can therefore be used as an effective two waveband pho- 
tometric redshift. This technique is important for current and fu- 
ture large area photo metric surveys which hope to study cluster 
abundance (e.g. RCS.lGladders & Yedl200ol: Pan- STARS S ; S DS S , 
Adelman-McCarthv et alj 120061 ; UKIDSS. lLawrence et all l200l 
Swinbank et alj|2007l) . 

In this work we will investigate the red sequence slope evolu- 
tion in observed optical (S(V — I) /SI) and near infrared (S(J — 
K)/SK) bands and rest frame optical (S(U — V)/8V) for a ho- 
mogeneous sample of X-ray selected galaxy clusters in the range 
2 ~0-l. This is the most comprehensive study of the red sequence 
slope undertaken thus far. We compare our findings with latest syn- 
th etic slopes calculate d from analysis of the semi-analytical model 
of Bower et al. (2006) based on the Millennium N-body simulation 
dSpringel et alj|2005f) . This model includes feedback from active 
galactic nuclei (AGN) which quenches star formation in massive 
halos to match the observed break in the luminosity function seen 
at bright magnitudes. Another important process in this model is 
'strangulation' which describes the stripping of a galaxy's hot gas 
halo as it falls into a cluster leading to a cessation in star formation 
in lower mass galaxies thus providing a mechanism for the build up 
of passive red galaxies in the cluster ( lLarsonet all 19801) . 
Lambda Cold Dark Matter (ACDM) cosmology (Q M =0.3, 
Qvac =0.7, Ho =70 km s _1 Mpc -1 ) is used throughout this work. 



2 OBSERVATIONS AND REDUCTION 
2.1 Cluster sample 

To observe the evolution of the red sequence we study X-ray se- 
lected clusters in the range z =0-1 in both near-infrared and optical 
bands. This is a large sample of the most X-ray luminous clusters 
known (Lx > 10 44 erg s" 1 , 0.1 - 2.4 keV) which, with the excep- 
tion of a handful of additional archival clusters, are all sourced from 
the ROentgen SATellite (ROSAT) All Sky Survey, although belong 



to the various sub-samples named below. Such massive clusters are 
ideal for this study as their colour-magnitude relations are well pop- 
ulated. The motivation for studying an X-ray selected sample of 
clusters is to ensure that we are observing objects in similar high 
mass, high density environments. This homogeneity is key to our 
study as we wish to compare clusters over a range of redshifts. 

The main cluster sample s studied in this pa p er, the MAs- 
sive Cluster Survey (MACS, lEbeling et ail 1200 ll ; lEbeling et ail 
2007h and the Las C ampanas/AAT Rich Cluster Survey (LARCS, 
Pimbblet et ^]|200ll . 2006), are sufficiently X-ray luminous that 
they should correspond to the most extreme environments at their 
respective epochs. The median X-ray luminosities of the high and 
low redshift samples are 16.1 and 7.4xl0 44 erg s _1 , respectively, 
corresponding to a difference of less than a factor of 2 in the typ- 
ical total mass dPopesso et alj|200^) . However, an important is- 
sue to address is that the mass of the z~0.5 progenitors of the 
LARCS clusters may be even lower than the MACS clusters. If 
we include the gr owth in halo mass through N-body mergers from 
Bowe r et al. I d2006h . we see that the progenitors of the LARCS clus- 
ters at z ~0.5 could be up to 3 times less massive than the MACS 
sample (with corresp onding X-ray luminosities of 2 x 10 44 erg s _1 , 
iPopesso et al.ll2005l) . There is no evidence for strong variations in 
the galaxy luminosity function between clusters spanning such a 
relatively modest difference in typical mass (de Propris et al. 1999). 
We therefore expect that any differences between the galaxy pop- 
ulations in these two samples will primarily reflect evolutionary 
differences between z ~0.5 and z ~0. 1. 

The 25 clusters observed in V and I or B and R bands belong 
to the following surveys: MACS (Vf555W and Vfsuw, Cycle 12 
GO project 9722); LARCS (B and R) and archival Hubble Space 
Telescope (HST) data (V FS5SW / F606W and Vf&uw)- The obser- 
vations were taken with the instruments: the Advanced Camera 
for Surveys (ACS) and the Wide Field Planetary Camera (WFPC) 
on HST and the lm Swope telescope at Las Campanas Observa- 
tory. For a more detailed descri ption of the da t a and reduction see 
IStott et al.l d2007l) (MACS) and lPimbblet et al.l i l200lh (LARCS). 

The 35 clusters studied in the near-infrared (J and K band) 
belong to: The MACS surve y ( lEbeling et ail 1200 ll) . the LARCS 
survey dPimbblet et all 1200 ll) . th e ROSAT Brightest Cluster Sur- 
vey (BCS) and extended BCS l Ebe ling et alTboOOh and the X - 
ray Brightest Abell Clusters Survey (XB ACS, lEbeling et aljl996t) . 
The observations were taken with: the Wide field InfraRed Cam- 
era (WIRC) instrument on the Palomar 200" Hale telescope the 
Infrared Spectrometer And Array Camera (ISAAC) on the Very 
Large Telescope (VLT) and the Wide Field CAMera (WFCAM) 
on United Kingdom Infrared Telescope (UKIRT). An additional 
low redshift data point is included for the Coma Cluster which we 
sourced from a combination of the Two Micron All Sky Survey 
(2MA SS) extended and poi nt source catalogu es dSkrutskie et al.l 
2006). We refer the reader to lStott et al] d2008l) for a more detailed 
description of the data and their reduction. 

We include a high redshift cluster (C1J1226.9+3332, z ~ 
0.9) from the W i de Angle ROSAT Pointed Surveys (WARPS, 
IScharf et al.ll 19971 : IJones et alj|l998h . This X-ray selected cluster 
was observed in the J and K ban ds using the U KIRT Fast-Track 
Imager (UFTI) camera on UKIRT dEllis & Jonesll2004h . 

2.2 Photometry 

The colour (V — I, B — R and J — K) photometry extracted 
for the cluster members employs 9kpc diameter apertures for 
colours and the magnitude used is SEXTRACTOR's 'Best' magni- 
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tude teertin & Arnoutslll996h . This choice of aperture is signifi- 
cantly greater than the seeing conditions in which the low redshift 
ground-based optical and near-infrared data were taken, typically 
~1.0 - 1.2 arcsec, giving LARCS sample photometric apertures of 
~4" whereas the higher redshift MACS sample is extracted with 
~1.4" apertures, considerably larger than the HST ACS PSF. We 
are therefore confident that we are collecting a comparable fraction 
of light for galaxies at different redshifts. We observe no trend be- 
tween colour and seeing for our LARCS low redshift optical data, 
although the seeing is consistent, and as the high redshift sample 
is observed with HST this is not a concern. We also investigate this 
for the z ~ 0.1 — 0.3 clusters observed in the near-infrared wave- 
bands for which colour is again independent of seeing, even though 
the conditions are more variable for our additional ROSAT clusters 
sample (0.9"- 1.5"). 

The SEXTRACTOR colour photometry was run in dual mode, 
on PSF matched images using the IRAF PSFMATCH package, with 
the 9kpc 'red' -band (R, I, K) apertures used to extract the cor- 
responding 'blue'-band (B, V, J) photometry. This is to ensure 
the same size aperture for both bands which is important for good 
colour determination. Star and Galaxy separation is performed us- 
ing SEXTRACTOR where detected objects with CLASS.STAR<0.1 
and/or J ~ K >0.95 (for the near-infrared sample) are classified 
as galaxies with the stars removed from the analysis. We only con- 
sider galaxies within 600kpc radius of the cluster centre, to limit 
contamination from field galaxies. All of our observations reach a 
depth which allows us to see to at least 4 magnitudes fainter than 
the BCG and thus perform a reliable fit to the slope and probe the 
sub-L* galaxy population where the red sequence build up is taking 
place. 

A potential problem with all aperture photometry is that the 
galaxies themselves may have significant internal colour gradients 
thought to be associated with a metallicity gradient between the 
outer and inner parts of the system. As stated above we are observ- 
ing the same fraction of light from galaxies at all redshifts which 
should limit this effect unless there is a significant evolution in the 
colour gradient of galaxies with redshift, although this would be an 
interesting effect in itself. Studies of the internal colour gradients of 
moderate redshift cluster galaxies using high quality //ST data have 
found no such evolution in mean colour gradient and redshift, with 
the indiv idual cluster galaxies a ppearing to have random colour 
gradients iTamura & Ohtal | |2000|) . We therefore conclude that any 
variations in internal colour gradient would act randomly to make 
galaxies appear either redder or bluer resulting in increased scatter 
of the red sequence but not an evolution in its slope. 

2.3 WFCAM data 

We obtained near-infrared J and K band data for 8 of the LARCS 
z ~ 0.1 clusters with WFCAM on UKIRT. WFCAM consists of 4 
detectors in a square, each separated by a gap comparable in size 
to a single detector, with a central autoguider. Each detector is a 
Rockwell Hawaii II 2048 x 2048 PACE HgCdTe array, with pixel 
size 0.4 arcsec. WFCAM single pointings or 'footprints' do not ob- 
serve a contiguous area of sky so to create a mosaiced image four of 
these WFCAM footprints are tiled to make a 4 detector x 4 detec- 
tor image. In addition to this tiling the WFCAM images incorporate 
'microstepping' which is a small dither pattern in the observation. 
The purpose of this process is to improve the point spread function 
sampling of the observations as the 0.4" pixel size of WFCAM is 
comparable to the best atmospheric seeing. The microstepping is 
performed by observing 4 images each offset from each other by 



a whole number of pixels and a 1/2 pixel in a 2 x 2 grid pattern. 
These images are then coadded together with the fractional offsets 
taken into account so that the resulting image has double the spatial 
sampling. Therefore, a detector image with an original resolution of 
0.4 "/pixel and size 2048 x 2048 pixels are converted to a higher 
resolution 0.2 "/pixel 4096 x 4096 image. 

The observations took place during service mode on the nights 
of 5th December 2006 and 19th of April 2007 in 1.0" and 1.2" see- 
ing conditions respectively. The total integration times for the J 
and K band images are 200s each, composed of 10s exposures in 
a 5 point dither pattern with 2x2 microstepping. To create a con- 
tiguous Mosaic of WFCAM images we run the TERAPIX SWarp 
software on the data which gives an image with a world coordi- 
nate system consistent to within 0. 1" of the 2MASS and the United 
States Naval Observatory (USNO) star catalogue across the whole 
field of view. To extract the J and K band catalogues from the mo- 
saiced images we use the SEXTRACTOR software in dual mode so 
that the K band catalogue detections are used to extract the J band 
photometry in an identical process to above. The photometry is cal- 
ibrated using the 2MASS point source catalogue and we find the J 
and K band 5 sigma vega limits of these data are 19.50 and 17.75 
magnitudes respectively. 

2.3.1 Quantifying the dwarf to giant ratio of the LARCS clusters 

The WFCAM instrument allows us to create contiguous images 
that are ~0.9 degrees on the side which corresponds to 5.6Mpc at 
2 = 0.1. With images of this size we can study both the cluster and 
the surrounding field galaxy population. A number of recent papers 
have studied the form of the colour-magnitude relation in galaxy 
clusters with some repo rting a dearth of faint red sequence galax- 
ies at high redshift (e . g.lKodama et al. 2004; De Lucia et alj|2004 
|Pe Lucia et"al]|2007l) . lstott et alj d2007F demonstrate that faint red 
galaxy population in clusters is built up over time by analysing the 
ratio of faint to luminous galaxies along the red sequence in a ho- 
mogenous sample of clusters at z ~ 0.1 and z ~ 0.5. One of the 
limiting factors in that study was the uncertainty in the faint end sta- 
tistical field correction for the LARCS z ~ 0.1 colour-magnitude 
relations. However the uncertainty in the field correction is reduced 
if near infrared observations are employed as it is easier to isolate 
the cluster red sequence. A comparison of Fig. Q] with the equiva- 
le nt optical colour m agnitude diagram, Fig. 1 (bottom right panel), 
of lStottetal] ( l2007l) illustrates this point. 

By ana lysing the WFCA M data in an identical way to that 
described in lStott et all d2007h in concert with the UKIDSS Deep 
eXtragalactic Survey (DXS, Survey Head: Alastair Edge) field ob- 
servations we can successfully quantify both the cluster and field 
galaxy population. The field correction is performed by dividing 
the colour-magnitude space of the cluster and field samples into a 
two-dimensional histogram a nd subtracting th e field number counts 
from the cluster population jStott et ai]|2007l) . The Poisson errors 
from this field correction are folded through into the final result. 
To determine the relative numbers of faint and luminous galax- 
ies on the red sequence we define the red sequence dwarf to gi- 
ant ratio (RDGR). This quantity is defined as the number of dwarf 
(faint) galaxies divided by the number of giant (luminous) galax- 
ies on the cluster red sequence after the statistical subtraction of 
the field galaxy population. The divid ing line betw een dwarfs and 
giants, for consistency with the IStott et al] d2007l) result, is given 
in absolute V band magnitude, Mv = —19.9 and a limiting mag- 
nitude of Mv ~ —17.75 for the faintest dwarfs, a correction for 
passive evolution is applied when studying high redshift samples. 
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Table 1. Details of the cluster sample used in our analysis w hich, with the excep tion of a handful of additional clusters, belong to ROSAT All Sky Survey. The 
LARCS and MACS redshifts are from lPimbblet et al.| (20061) and lEbeling et alj {2007.} respectively. Kopt. 

and kjk are the observed optical and near-infrared 

slopes. 



Cluster R.A. Dec. z Lx Kapt. Kjk 

(J2000) (lO^ergs" 1 ) 



LARCS z ~ 0.1 Sample 



Abell 22 




00 20 38.64 


-25 43 19 


0.142 


5.3 


-0.072 ±0.014 


-0.023±0.010 


Abell 550 




05 52 51.84 


-21 03 54 


0.099 


7.1 


-0.055 ±0.013 


-0.024±0.011 


Abell 1084 




10 44 30.72 


-07 05 02 


0.132 


7.4 


-0.058 ±0.013 


-0.033±0.007 


Abell 1285 




1 1 30 20.64 


-14 34 30 


0.106 


5.45 


-0.036 ±0.020 




Abell 1437 




12 00 25.44 


+03 21 04 


0.134 


7.7 


-0.037 ±0.017 


-0.010±0.009 


Abell 1650 




12 58 41.76 


-01 45 22 


0.084 


7.8 


-0.047 ±0.01 1 


-0.019±0.009 


Abell 1651 




12 59 24.00 


-04 11 20 


0.085 


8.3 


-0.048 ±0.009 


-0.002±0.007 


Abell 1664 
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Abell 2055 




15 18 41.28 


+06 12 40 


0.102 


4.8 


-0.066 ±0.01 1 


-0.021±0.009 


Abell 3888 




22 34 32.88 


-37 43 59 


0.153 


14.5 


-0.047 ±0.022 




MACS z = 0.4-0.7 Sampk 














MACS J0018 


5+1626 


00:18:33.68 


+ 16:26:15 


0.541 


18.74 




-0.048±0.009 


MACS J0025 


4-1222 


00 25 15.84 


-12 19 44 


0.478 


12.4 


-0.106 ±0.021 


-0.047±0.010 


MACS J0257 


6-2209 


02 57 07.96 


-23 26 08 


0.504 


15.4 


-0.105 ±0.014 


-0.065±0.012 


MACS J0454 


1-0300 


04:54:11.13 


-03:00:53.8 


0.550 


16.86 




-0.040 ±0.013 


MACS J0647 


7+7015 


06 47 51.45 


+70 15 04 


0.584 


21.7 


-0.090 ±0.024 




MACS J0712 


3+5931 


07 12 20.45 


+59 32 20 


0.328 


6.8 


-0.044 ±0.006 




MACS J0717 


5+3745 


07 17 31.83 


+37 45 05 


0.548 


27.4 


-0.068 ±0.012 




MACS J0744 


8+3927 


07 44 51.98 


+39 27 35 


0.686 


25.9 


-0.081 ±0.023 




MACS J0911 


2+1746 


09 11 10.23 


+ 17 46 38 


0.506 


13.2 


-0.087 ±0.011 




MACS Jl 149 


5+2223 


1149 34.81 


+22 24 13 


0.544 


17.3 


-0.101 ±0.011 




MACS J 1354 


6+7715 


13 54 19.71 


+77 15 26 


0.397 


8.2 


-0.057 ±0.018 




MACS J1359 


8+6231 


13:59:54.32 


+62:30:36.3 


0.330 


8.83 




-0.039 ±0.007 


MACS J 1423 


8+2404 


14 23 47.95 


+24 04 59 


0.544 


15.0 


-0.078 ±0.013 




MACS J2129 


4-0741 


21 29 25.38 


-07 41 26 


0.570 


16.4 


-0.115 ±0.017 


-0.068±0.009 


MACS J2214 


9-1359 


22 14 56.51 


i a nn 1 1 

— 1 <+ UU 1 / 
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Additional ROSAT All Sky Survey clusters 












Abell 1 15 




00:56:00.24 


+26:20:31.7 


0.197 


14.59 




-0.031 ±0.006 


Abell 209 




01:31:52.51 


-13:36:41.0 


0.209 


13.75 




-0.022 ±0.009 


Abell 291 




02:01:46.80 


-02:11:56.9 


0.196 


4.24 




-0.024 ±0.009 


Abell 665 




08:30:57.34 


+65:50:31.4 


0.182 


16.33 




-0.020 ±0.008 


Abell 773 




09:17:53.57 
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Abell 1201 




11:12:54.50 


+ 13:26:08.9 


0.169 


6.28 




-0.013 ±0.007 


Abell 1246 




11:23:58.75 


+21:28:47.3 


0.190 


7.62 




-0.031 ±0.008 


Abell 1703 




13 15 00.70 


+51 49 10 


0.258 


8.7 


-0.057 ±0.012 




Abell 1758 




13:32:38.59 


+50:33:38.7 


0.279 


11.68 




-0.027 ±0.01 1 


Abell 1763 




13:35:20.14 


+41:00:03.8 


0.223 


14.93 




-0.029 ±0.010 


Abell 1914 




14:25:56.64 


+37:48:59.4 


0.171 


18.39 




-0.021 ±0.006 


Abell 21 11 




15:39:41.81 


+34:24:43.3 


0.229 


10.94 




-0.033 ±0.012 


Abell 2163 




16:15:33.57 


-06:09:16.8 


0.203 


37.50 




-0.022 ±0.006 


Abell 2218 




16:35:49.39 


+66:12:45.1 


0.176 


9.30 




-0.025 ±0.007 


Abell 2445 




22:26:55.80 


+25:50:09.4 


0.165 


4.00 




-0.027 ±0.009 


RXJ1720.1 + 


2638 


17:20:10.08 


+26:37:33.5 


0.164 


6.66 




-0.023 ±0.007 


Zw 1432 




07:51:25.15 


+ 17:30:51.8 


0.186 


5.27 




-0.022 ±0.010 


Additional archival clusters 














C1J0152-1357 


01 52 43.91 


-13 57 21 


0.831 


5.0 


-0.084 ±0.026 




C1J1226.9+3332 


12 26 58.13 


+33 32 49 


0.890 


20.0 


-0.088 ±0.016 


-0.088±0.010 


Coma Cluster 


12:59:48.70 


+27:58:50.0 


0.0231 


7.26 


-0.075 ±0.015 


-0.017 ±0.009 


MS1054-0321 


10:57:00.20 


-03:37:27.4 


0.830 


23.30 




-0.073 ±0.01 1 


RCS0224-0002 


02:24:00.00 


-0:02:00.0 


0.770 


0.70 




-0.055 ±0.020 
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Figure 1. The (J — K) vs K colour-magnitude diagram for the cluster 
Abell 1084 which demonstrates the ease at which near-infrared colours sep- 
arate the red sequence from stellar and field galaxy contamination 
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Figure 2. The evolution of the of the red sequenc e giant dwarf ratio (RDGR) 
with redshift comparing the MACS sample from lStott et"aT] fe007l) with the 
LARCS WFCAM sample. A fit of the form (1 + z)-P to the LARCS and 
MACS samples is plotted which yields (3 = —2.1 ± 0.3. 



We convert these limits to apparent K band magnitudes, K ~ 14.7 
and K ~ 16.8 respectively (de pending on the precise redshift), us- 
ing a lBruzua l & Chariot] j2003l> simple stellar population with solar 
metallicity and a formation redshift, Zf = 5. 

The weighted mean of the RDGRs of the 8 LARCS clusters 
from this sample is 2.59±0.24 compared to the lStott et ail J2007I) 
result of 2.93±0.45. We can now say that when compared to our 
z ~0.5 MACS sample, where the DGR=1.33±0.06, the ratio of 
dwarf to giant galaxies on the cluster red sequence has increased 
by a factor of 1.95±0.20 in the past 5 Gyr which confir ms and im- 
proves the significance of the result o f lStottetal.l {2007) from 3a to 
5a. This is strong evidence for a significant build up of the red se- 
q uence in ma s sive c lusters. Fig.|2]is the equivalent of Fig. 2 (right) 
of IStott et aT] d2007h and displays the evolution in the RDGR with 
redshift for rich galaxy clusters. To parameterise this evolution we 
fit a (1 + z) _/3 power law to the LARCS and MACS samples where 
/3 = — 2.1 ±0.3 with all clusters consistent with the fit. 



3 ANALYSIS AND RESULTS 
3.1 Fitting red sequence the slope 

Figures [3] and [4] are examples of the prominent red sequences we 
see in the rich clusters of our sample. To fit the slope of the red se- 
quence we need to use a robust and consiste nt method. We therefor e 
employ the same technique as described in Glad ders et al.l dl998l) . 
an iterated 3 a clipped fit. The fit is performed as follows: first we 
set a limiting magnitude for the red sequence which corresponds 
to the mean of the next 2 brightest galaxies down from the BCG 
magnitude+3. The reason for this is that in some cases there can 
be a large luminosity gap (~1 mag) between the BCG and the start 
of the red sequence so this ensures we measure the slope without 
this gap. The fit is performed over 3 magnitudes to incorporate the 
sub-L* galaxy populati on which show t he strongest evidence for 
red sequence build up JStott et al]|2007l) . We select the region of 
colour magnitude space containing the red sequence and estimate 
an initial fit from visual inspection. The residuals about this esti- 
mate are calculated and a Gaussian is fitted to the resulting colour 
distribution with the slope removed. The peak of this Gaussian cor- 
responds to the red sequence. We then perform a fit to points that 
are within 3 sigma of this fit. This is a two parameter linear fit of 
the form y = kx + c where k is the slope of the red sequence. 
The process is iterated until it con verges to a solution. We confirm 
the work of Gladders et al. d 19981) that this is a robust method for 
reasonable choices of limiting magnitude and initial fit. 

It should be noted that our colour-magnitude diagrams are not 
field corrected as the rich clusters in our sample have well pop- 
ulated red sequences in contrast to the field. At the colours and 
magnitudes we consider this a very small contribution to the red 
sequence, typically 5% contamination for the MACS clusters and 
10% for the LARCS clusters. We perform a statistical field correc- 
tion test to a sub-sample of our low z clusters and find that the slope 
of the sequence varies randomly by less than la than that obtained 
for the uncorrected sequence. We therefore feel justified in not ap- 
plying this correct ion. For details of the statistical field correction 
technique used see IStott et alj J2007h . 



3.2 Model slope evolution 

We form a prediction of how the red s equence slope evol ves, by 
analysing the semi-analytical model of Bower et al. ( 2006) based 
on the Millennium N-body simulation JSpringel et alj|2005h . We 
access these data online through Virgo Millennium Database. The 
Bowe r et al. liEooi) model accounts for the recent observations that 
the stellar mass in bright galaxies was in place at high redshift by 
including feedback from active galactic nuclei. It also includes a 
prescription for the build of the red sequence via 'strangulation'. 
This describes a process whereby galaxies falling into a cluster are 
stripped of their hot gas reservoir upon int eraction with the intra- 
cluster medium quenching star formation ^Larson et all 198(3) . 

We define model clusters as dark matter haloes above a thresh- 
old mass which contain bound sub-haloes (galaxies) with obser- 
vational properties assigned to them by semi-analytic modelling. 
To allow comparison between models and observation we as- 
sume that the model cluster members are analogous to the red se- 
quences of our obs erved clusters. We need to select clusters from 
iBower et all J2006l) which are comparable in mass to those we ob- 
serve. By using the relation between cluster mass and X-ray lumi- 
nosity (Po pesso et alj |2005) we include only the simulated haloes 
with M200 > 3.4 x 10 14 M Q which corresponds to Lx > 10 44 erg 
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Figure 3. The LARCS optical colour-magnitude diagrams with the red sequence fits overplotted. Representative error bars are shown for a range of magnitudes. 
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Figure 4. The main MACS sample optical colour-magnitude diagrams with the red sequence tits overplotted. Representative error bars are shown for a range 
of magnitudes. 
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s _1 . To select only the most massive systems for comparison with 
our observations we limit our red sequence fit to a stack of the top 
five ranked clusters by mass in each redshift bin. The model slope 
is robust to our choice of lower mass limit as we see no signifi- 
cant change in our slope values when using all haloes with masses 
greater than mean halo mass in the required redshift interval. 

We model the red sequence slope ev olution by creating 
stacked colour magnitude diagrams from the iBower et alj J2006T) 
model output at a distinct set of redshift intervals between z = 
and z — 1. For the creation of these colour-magnitude diagrams we 
ensure we only study the passive red sequence galaxies by selecting 
galaxies with no current star formation (L(Ha) =0). The models 
provide no spatial information but we assume that the synthetic 
colours we calculate from the total magnitudes can be compared 
with the observed aperture magnitudes allowing for a normalisation 
between the model and observed slope evolution at low redshift. 

We calculate the model for the observed frame slope evolution 
and its errors by fitting to the stacked synthetic red sequence slopes 
at each redshift interval with the method described in As with 
the observations, the model slope is shown to steepen with redshift. 

3.3 Observed slope evolution 

Fig. [5] displays the observed red sequence slope evolution for our 
near infrared sample. The slope (kjk = S(J — K)/8K) is shown 
to steepen with redshift. This steepening will have contributions 
from both K correction and perhaps an age or a metallicity evolu- 
tion. The contribution from K correction is due to the observed J 
and K bands sampling increasingly bluer rest wavebands at higher 
redshift which affects galaxies differentially along t he red sequence 
depending on their spectr al energy distributions ( Gladde rs et al.l 
1 19981: iPimbblet et af]|200lh . 

When the model described in <j3 .21 is plotted with our com- 
plete near-infrared dataset in Fig[5]we find good quantitative agree- 
ment between the two with the data having an rms scatter about the 
model of 0.009 which is comparable to the calculated error on the 
model (mean la error is 0.006). 

We now investigate the slope (kvi = 8{V — I) /SI) evo- 
lution for our optical observations. As above, the simulated slope 
evolution plotted is calculated from semi-analytical model of 
lBoweretalJ ( l2006l) . There is an additional complication as our op- 
tical data are sourced from two different filter sets so we have to 
account for the observed difference in red sequence slope between 
them. To achieve this we normalise all data to the 8(V — I)/ SI 
filters by correcting for the difference between the S(V — I) /SI 
m odel and the model fo r the S(B — R)/8R filter combination (as 
in iGladders et ai¥ l998). The resultant data-points and model are 
plotted in Fig. [6] The slope is shown to increase with redshift as 
above. The rms scatter about the model is 0.017 (compared to the 
mean 1 a error in the model of 0.008) so as with the near-infrared 
observations we find good agreement between the model and the 
data. 



3.3.1 Rest frame slope evolution 

To investigate the intrinsic evolution of the red sequence slope we 
need to study clusters observed with matched rest frame photom- 
etry. This is to quantify the proposed contribution from the build 
up of the red sequence without the additional effect of K correc- 
tion. We can observe this intrinsic slope evolution as our two main 
sets optical observations, the z ~0. 1 LARCS clusters in (B — R) 



colour and the z ~0.5 MACS clusters in (V — J) colour, both 
correspond to the rest frame (U — V) colour at their respective 
epochs. (U — V) colour straddles the 4000A break and therefore 
is a good discriminant between the red sequence and star forming 
cluster members/foreground galaxies. We confirm that the filters 
are well matched to rest frame (17 — V) as we find a colour term of 
1.0 ±0.05 between the B - R at z =0.11 and V - I at z =0.53 
using the technique described in Bla keslee et al.l (2006). 

The evolution of this intrinsic slope (kuv = S(U — V)/SV) 
is plotted in Fig. [7j We include an additional low redshift data- 
point for the C oma cluster calculated from the Sloan Digital Sky 
Survey (SDSS. lAdelman-McCarfhv et alj|2006l) u and g filter pho- 
tometry. In this figure we can see that the intrinsic optical slope, 
Kuv, does evolve with redshift as the intermediate z MACS clus- 
ters have a steeper red sequence than their low z LARCS counter- 
parts. The weighted mean values for the LARCS and MACS Kuv 
are -0.053±0.004 (s.e.m) and -0.092±0.004 (s.e.m) respectively, a 
difference of 6.5 o. We can therefore say that there is a real con- 
tribution to the slope evolution from factors other than K correc- 
tion. The fit to the data in Fig. [7] is of the form (1 + z) 13 where 
(3 — 1.77 ± 0.25. We note that the Coma Cluster has a steeper 
slope than the rest of our low redshift sample, although only a ~2cr 
discrepancy from the fit, which we may expect as it is found to 
have lower than average dwarf-to-giant ratio along its red sequence 
sugges ting it is still und ergoing faint end and therefore slope evo- 
lution ( Stott et al. 2007). Including the Coma Cluster does not have 
a significant affect on the fit, (3 — 1.67 ± 0.26. 

The dashed line plotted on Fig. UJ is the rest frame kuv = 
S(U — V)/ 8 V slo pe calculated from the semi-analytic model of 
Bowe r et al. I d2006h . normalised to the mean value of the LARCS 
sample slope, which shows only a mild evolution with redshift and 
is therefore unable to replicate the rest frame slope change ob- 
served. From this we can conclude that the major contribution to 
the agreement between the observed slope evolution in i j3.3l and 
that derived from the semi-analytic modelling shown in Figs[5]&[6] 
was the K correction differential between high and low mass galax- 
ies along the synthetic red sequence and that a significant, intrinsic, 
slope evolution is not predicted by the models. 

3.4 Evolution with other observables 

We now investigate whether there are further trends in red sequence 
slope with other observable cluster properties to ensure that the re- 
sult seen in ^3.3. 1 l is not due to a secondary correlation. The most 
obvious of these being the X-ray luminosity, which is a proxy for 
mass in a relaxed system. In Figs [8] and [9] we plot the residual val- 
ues of the slope about the model line in Fig. [5] and fit in Fig. UJ 
against Lx ■ From this we see no significant trend between scatter 
about the model and Lx as the Pearson correlation coefficients, r 
, for the near-infrared and optical data are 0.17 and 0.24 respec- 
tively which is actually a weak anti-correlation between steepening 
negative red sequence slope and Lx - However, because of the mag- 
nitude of the errors on both the individual slope measurements and 
the models, we need to quantify what level of trend the errors could 
accommodate before becoming observable. For the near-infrared 
sample we find that a steepening greater than 0.010 in the slope in 
Lx = 26.6 x 10 44 erg s _1 can be ruled out with a 3 sigma con- 
fidence whereas for the rest frame (U — V) sample we find this 
magnitude of change can occur in Lx = 12.3 x 10 44 erg a -1 . The 
difference in mean Lx between the LARCS and MACS samples 
is 8.5 x 10 44 erg s _1 for which we can rule out, with a 3 sigma 
confidence, increases in slope greater than 0.007, if we use the 
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Figure 5. The evolution of the red sequence slope (re jx ) for our near-infrared sample. A model calculated from Bower et al. ( 2006) is included for comparison 
with theory 




(U — V) result, or 0.003 if we use the J — K result. For the rest 
frame 8(U — V)/SV slope evolution in Fig. [7] we see that the ob- 
served change in rest frame slope of ~ 0.04 cannot be accounted 
for by the difference in Lx alone and is, we believe, a definite 
trend with redshift. In Fig. [10] we also show that there is no sig- 
nificant trend between red sequence slope and another mass proxy, 
the cluster velocity dispersion, a (r= 0.27), which is again a weak 
anti-correlation. As above we can rule out increases in slope greater 
than 0.010 with 3 sigma confidence between the average velocity 
dispersions of the MACS and LARCS samples. 



In addition to the X-ray luminosity and a, we can look to the 
degree of BCG dominance as an indicator of the local environment 
within the cluster core. This parameterises the luminosity gap be- 
tween the BCG and the next brightest galaxies on the red sequence 
and is defined as Ami_2,3 = (rri2 + mz)j1 — m\ where m\ is 
the magnitude of the BCG and m.2 and mz are t he magnitudes of 
the 2nd and 3rd brightest members respectively 1 Ki m et al]|2002l : 
IStott et al.l 12008), The BCG may be the dominant elliptical in a 
cluster centre containing much smaller galaxies or it may be in a 
system were it is only marginally brighter than the next brightest 
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Figure 7. The evolution of the rest frame red sequence slope (kjjy) for our optical sample. The triangle is the data point for the Coma Cluster. The solid line 
is a fit to the MACS and LA RCS data of t he form (1 + z)" where /3 = 1.77 ± 0.25. The dashed line is the rest frame slope evolution calculated from the 
semi-analytic model of lBoweretaT]j2006l) . 
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Figure 8. The residuals about the models in Fig.[5]plotted against X-ray lu- 
minosity. Filled circles and squares represent galaxy clusters with redshifts, 
2 < 0.2 and z > 0.2 respectively, to demonstrate there is no redshift 
dependence. The correlation coefficient, r, for this plot is 0.17. 



members. In Fig. [TT]we demonstrate that there is no strong trend 
between red sequence slope and the degree of BCG dominance 
(r= —0.09) and can rule out increases of slope greater than 0.01 
in 0.61 magnitudes of dominance with a 3 sigma confidence. 

The above results suggest that different cluster environments 
do not have a strong effect on the near-infrared or optical red se- 
quence slope at a given redshift, demonstrating that the result in 
^13.31 is robust. Previous studies within a similar Lx range have 
also seen homogeneity in other related cluster properties such as 
the shape of the luminosity function an d the blue galaxy fraction 
jde Propris et alj|l999l ; lwake et al.ll2005b . 
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Figure 9. The residuals about the models in Fig.|7]plotted against X-ray lu- 
minosity. Filled circles and squares represent galaxy clusters with redshifts, 
z < 0.2 and z > 0.2 respectively, to demonstrate there is no redshift 
dependence. The correlation coefficient, r, for this plot is 0.24. 

4 DISCUSSION 

In this work we have found a significant evolution in the rest frame 
slope of the red sequence in rich galaxy clusters between z ~ 0.5 
and z ~ 0.1. We propose that this intrinsic evolution is due to 
galaxies f alling into the cluster core and transfor ming onto the red 
sequence dDe Lucia et al. 2007; Stott et al.ll2007l) . If th ese galaxies 
have undergone recent star formation in filaments (e.g. lPorter et al.l 
2008), which has been quenched by interactions as they fall into 
the cluster, they will appear bluer than other passive galaxies in the 
cluster due to their young age. These galaxies will redden more 
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Figure 10. The residuals about the models in Fig. [7] plotted against veloc- 
ity dispersion, a. Filled circles and squares represent galaxy clusters with 
redshifts, z < 0.2 and z > 0.2 respectively, to demonstrate there is no red- 
shift depende nce. The MACS and LA RCS cluster velocity d ispersions are 
sourced from Ebeling et al. (2007) and Pimbblet et al. (2006) respectively. 
The correlation coefficient, r, for this plot is 0.27. 
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Figure 11. The residuals about the models in Fig. \5\ plotted against BCG 
dominance. Filled circles and squares represent galaxy clusters with red- 
shifts, z < 0.2 and z > 0.2 respectively, to demonstrate there is no redshift 
dependence. The correlation coefficient, r, for this plot is -0.09. 
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dl998h . iBlakeslee etafl d2003r) and iMei et all j2006t) who find no 
evidence for intrinsic slope evolution. There are several potential 
reasons for this disagreement. We first note the scatter in the values 
of the slope which could obscure the intrinsic evolution in previ- 
ous studies without the coverage per redshift bin we have in our 
sample. This scatter is important when choosing a low redshift 
comparison cluster, some earlier works use a single cluster, usu- 
ally Coma, which we have shown to have an unusually steep red 
sequence compared to the mean of our low z sample. The use of 
Coma would therefore tend to bias the result towards no evolu- 
tion. The choice of (U — B) rest frame colour in previous studies 
means that they are observing a smaller dynamic range in colour 
than our (U — V) sample. This is particularly important if the rest 
frame 'U — B' colours do not e ffectively straddle the 4000A break 
(e.g. figure 4 o flMei et"al . 2006) which may result in contamination 
of the red sequence by both foreground galaxies and star forming 
cluster members. We also ensure that we perform all of the red se- 
quence slope fits using the same homogenous method rather than 
importing slopes from other studies which may have a systematic 
offset. 

The good agreement found between our data and s lope evolu- 
tion m odels calculated from semi-analytical model of Bo wer et"al] 
d2006l) in Figures [5] and [6] is mainly due to the treatment of K 
correction rather than any significant slope evolution. This demon- 
strates that although such models include 'strangulation' of star- 
forming galaxies falling into cluster environments, they are unable 
to effectively reproduce the intrinsic evolution seen in our sample. 
This defic iency is also seen when studying the build up of the red 
sequence dGilbank & Baloghl2008l) . We speculate that an improved 
prescription for a strangulation-like process in the semi-analytic 
models will be an important factor in recreating observations of 
an intrinsic slope evolution and build up of the red sequence. 

When looking at slope trends with other observables we see 
no relationship between cluster slope and X-ray luminosity, ve- 
locity dispersion or BCG degree of dominance. This suggests that 
there is very little variation between the red sequence slopes due to 
the different cluster environments demonstrating that the main re- 
sult of this paper is robust to the differing global properties within 
our sample. This implies that searching for massive clusters us- 
ing the colour magnitude relation (CMR), as employed by cur- 
rent and future large area optical/near-infra red surveys, is a viable 
method (e.g. RCS, iGladders & Yeej l2000l : Pan-STARSS; SPSS 
Adelman-McCarthv et al.l 120061 ; UKIDSS. lLawrence etaf] 120071 
Swinbank et alj|2007h . 



rapidly with time than their old, luminous counterparts flattening 
the red sequence at progressively lower redshift. This relates to the 
concept of downsizing which describes the observation that star 
formation ceases in the largest galaxies fir st so that it takes p lace 
mainly in lower mass galaxies at late times dCowie et al .11 1996b . 

An alternative or additional explanation is that this result is 
caused by a differential chemical evolution along the red sequence. 
A differential chemical evolution along the sequence which pref- 
erentially enriches the faint galaxies could lead to a steeper red 
sequence slope at high redshift. The supposition of a metallic- 
ity differential along the sequence that decreases with redshift is 
borne out in the semi-analytic models as the gradient at z =0 
is -0.0029±0.0004 Z mag -1 and at z=l it is -0.061 ±0.009 Z 
mag -1 . However, this may be non-physical with the primary goal 
of semi-analytic model is to match the observed properties of the 
local Universe and not necessarily the evolution from high redshift. 
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